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In the past many papers have appeared which simulated 
surface growth with different growth models. The results 
showed that, if models differed only slightly in their 'growth' 
rules, the resulting surfaces may belong to different univer- 
sality classes, i.e. they are described by different differential 
equations. In the present paper we describe a mapping of 
"growth rules" to differential operators and give plausibility 
arguments for this mapping. We illustrate the validity of our 
theory by applying it to published results. 



I. INTRODUCTION 

During the past few years the kinetic roughening of 
surfaces has become a field of increasing interest. In 
particular, many papers have appeared concerned with 
computer simulations of surface growth [see e.g. Rcf. |. 

In general, the surface is characterized by a height h 
appropriate to a c? dimensional substrate of size L. The 
width of the surface w{t, L) at a time t is characterized 



dh{x,t) 
dt 



F + ri{x,t) 



(2) 



by w(t, L) — Y [hY — (h) where the bar denotes an av- 
erage. If the resulting surface is self-afRne it can be rep- 
resented by a dynamical scaling law 



(1) 



where the function f{x) constant for a; — > oo and 
f{x) ~ x^ , with P — a/z as a; — > 0. The unit of time 
corresponds to depositing L particles. For models which 
do not contain vacancies or overhangs, so-called solid on 
solid models (SOS), this means that the average height 
h and the time t are identical. The exponents a, [3 and 
z determine which universality class the given model be- 
longs to. In the present paper we wish to examine the 
implications of some of the standard assumptions in the 
theory, and in particular to provide a mapping of pre- 
scribed rules of growth to the corresponding differential 
operators appearing in the associated stochastic growth 
equations. 



II. THEORY AND DISCUSSION 

A. Simple growth systems 

The dynamical evolution of a surface prior to any 
movement of the deposited particles is presumed to be 
described by the equation 



where is a noise term with zero mean (i.e. ^ = 0) and 
F is the flux rate of incident particles. For the particu- 
lar case, usually considered in the literature, that is a 
constant we obtain h = Ft. The last result is only true 
if there arc no vacancies in the system. A similar result 
(but with F replaced by a larger entity, F' say) would 
also hold, if the vacancy concentration stays constant 
over time, i.e. in this case we would again have h — F't. 
However, for such a situation the starting equation would 
not be equal to Eq.(||), and F' ^ F, the incident flux. In 
general, the form of the starting equation for situations 
where vacancies and overhangs occur is far from obvi- 
ous and is a problem we will return to later. Before we 
can answer such questions it is pertinent to address a far 
simpler problem. Namely, what type of deposition pro- 
cess is Eq. (H) applicable to anyway? One class of such 
processes is the so called SOS models in which particles 
are deposited randomly at lattice sites, i.e. a number is 
chosen at random and the number of particles at the site 
characterized by that number is increased by one. This 
is usually referred to as (pure) random deposition (RD) 
and corresponds to a constant flux F with a random noise 
term rj. (We refer to it as pure since the deposited atoms 
are not permitted to move.) It is well known that the 
interface width for (pure)RD increases indefinitely with 
time, i.e. the associated surface is not self-affine with_the 
consequence that the width itself does not saturate.EJ 

A variant of (pure) RD is to allow the particles to move 
after they have been deposited. Two types of movement 
are possible. One is 'horizontal' movement in which the 
height of the moving particle does not change, and the 
other is vertical movement in which the height of the 
particle does change. With regard to the latter there are 
two possible types. One is upward vertical movement in 
which the height of the particle increases and the other 
is downward vertical movement in which the height of 
the moving particle decreases. In order to complete the 
rules of growth with regard to the 'allowed' moves we 
need to specify the conditions under which horizontal or 
vertical movement is terminated. We will refer to these 
collectively as 'sticking rules'. 

In what follows we will show that what superficially ap- 
pear to be 'trivial' or 'obvious' rules of movement (often 
designed for convenience of implementation on a com- 
puter) can have extremely subtle implications for the re- 
sulting mathematics (and, in particular, the associated 
differential equations purporting to model the resulting 
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surface structure). Furthermore, we will show that the 
implications for the mathematics (even for the same set 
of rules) can be different dependent on the dimensionality 
of the surface dl under consideration. 

To be definite, and also for reasons of simplicity, we 
will focus attention initially on a d' = 1 surface i.e. a 
two dimensional structure (c? — 2) having a y-axis (the 
height) and an cc-axis (the coordinate defining the posi- 
tion along the surface). Let us consider first a variant on 
the pure RD model in which we allow both horizontal and 
vertical motion (such motion being described under the 
general heading of surface relaxation processes). Within 
the constraint of SOS, the allowed moves could be the 
following. Permit the particle to be displaced horizon- 
tally up to a maximum of / lattice sites provided always 
that there is another particle immediately below it. If 
such a situation always prevails, then at the end of the 
I moves leave the particle where it is. If, on the other 
hand, in the course of making these moves the particle 
encounters a gap, then allow it to move vertically down 
until it reaches a position where there is again a particle 
immediately below it at which point the particle move- 
ment is terminated. Similarly, if in the course of making 
up to / moves, it meets another particle on the same level 
all movement is terminated. 

It is clear that surface relaxation processes involving 
such SOS movements have the effect of changing the noise 
term 77 in Eq (^ to a new value rj' . Thus this situation 
could be described by the equation 



(3) 



Such an equation, although correct, is of little use from 
an analytical viewpoint because we do not know the form 
of the noise term r/. However the essential characteristic 
of SOS movement that we can make use of is that it 
leaves the average height of the surface unchanged. In 
effect one lattice site loses a particle and another lattice 
site (which could be anywhere up to a distance I away 
horizontally) gains one. This suggests that we write the 
noise term ?]' in the form 



dG 
dx 



(4) 



where G is any (well behaved) function we care to choose. 
The reason being that the act of averaging over a finite 
(but large) number of discrete particles is deemed as be- 
ing equivalent (in a continuum description) to the act of 
integrating and then dividing by the 'length' of the re- 
gion of integration. Consequently, expressing the change 
in the noise as a divergence means (at least in the limit 
that the length of the system becomes infinite) that the 
average of 77' is zero (just like 77). This then gives the 
formal result 

{h) ^ Ft 

which is known to be correct for SOS models. (One might 
note that, strictly speaking, this result is not 'exact' for 



a finite system since the divergence term integrates to 
a small but finite entity. Such 'errors' are an inherent 
feature of the modelling of finite systems by continuum 
differential equations). Thus, for SOS models invoking 
'surface relaxation' we expect the generalisation of the 
equation of motion Eq (||), for (pure) RD to be 



dh_^ dG 
dt dx 



V 



(5) 



where G is, at this stage, some unspecified function. A 
computer model of dL = 1 RD plus surface relaxation was 
evaluated by FamilyEI many years ago. According to this 
author the resulting surface was self-affine, i.e. was de- 
scribed by Eq (|l|) with the value of the parameters a, (3 
and z being independent of the maximum number of lat- 
tice sites I over which a particle was allowed to move 
horizontally l& point we will return to later). Further- 
more Familya found that the scaling properties of the 
surface (or equivalently, the values of a, (3 and z ) were 
consistent with the choice 



G^u 



dh 
dx 



(6) 



with V constant. 

The resulting equation 



dh_^ ^d^h 
dt dx^ 



(7) 



is well known in the literature and is referred to as the Ed- 
wards Wilkinson EW equation.^ It is important to note 
exactly how the mapping of the horizontal and vertical 
motion onto a differential operator occurs in Eq (|^). The 

entity I'^jr (with v positive) corresponds (as shown in 
Fig. 1^) to a movement of particles from the top of a 'hill' 
to the bottom of it, i.e a flattening of the 'hill'. In other 
words it is consistent with movements both horizontally 
and vertically downwards. This manifests itself clearly in 
that model of Familyo in which all possible downwards 
moves of a particle are allowed to occur, and the parti- 
cle moves to the nearest neighbour position only in the 
horizontal movement (i.e. I — 1). The form of the as- 
sociated surface is shown in figure ^. Examination of 
the latter shows that the 'vertical' separation of nearest 
neighbour columns is also small for this situation. Hence 
the model leads self-consistently (as a result of the rules 
of movement) to a situation where both small horizon- 
tal and downward vertical motion is the norm. Thus the 
situation were all possible downward motion is permit- 
ted to occur with 100 % probability is replicated in the 
mathematics by the differential form 7^§-4. 

This having b^en established we now return to the 
claim by Familytl that the parameters a, (3 and z are 
independent of the magnitude of L To see that this can- 
not possibly be true let us take a d! — \ surface of length 
L and allow the number of horizontal moves to be up to 
a maximum of L. For such a situation it is clear that we 
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will obtain layer by layer growth - the reason being the 
following. During the first monolayer coverage 'islands' 
will develop of size Li say. With increasing coverage it be- 
comes increasingly probable that the next particle added 
will be on one of these islands. However, if the particle is 
allowed L horizontal moves, it will migrate readily to the 
edge of the island, go over the side and then adhere to 
the rim of it, thus extending the size of the island. There 
are several important results to be deduced from these 
considerations. First, if we have a surface of total length 
L, then as long as we allow horizontal moves up to some 
maximum I say, we may well find a series of surfaces be- 
longing to the same universality class. However this can 
only be true as long as / is less than some fraction of the 
total 'length' L of the surface. Once / exceeds this frac- 
tion, an increasing number of islands will coalesce and we 
will very quickly obtain layer by layer growth. In short, 
dependent on the size of the system, universality classes 
will only be obtained if we permit a limited number of 
horizontal movements only. The precise number is ill de- 
fined at this stage, but typically in the literature this 
number is chosen to be one or two only. 

B. Surface Relaxation Involving Vertical Motion 
Downwards Only 

We consider next vertical motion. To this end we will 
examine first what is usually referred to in the literature 
as (pure) ballistic deposition (BD). In such a model one 
imagines the particles being deposited vertically onto a 
substrate in a random sequence. The rule of movement is 
to move vertically down until the particle meets another 
particle. The latter can be either directly underneath or 
at the side of the moving particle. In either event move- 
ment is then terminated. As is well known, such a model 
does have overhangs and vacancies present in it. Let us 
pretend for the moment that we do not know the equa- 
tion of motion for the evolution of such a surface and 
consider instead a variant on this model. Namely (pure) 
BD plus surface relaxation. The specification of the latter 
is the following. If there is another particle immediately 
beneath the one under consideration then no movement 
occurs. On the other hand, if there is no particle im- 
mediately beneath the one under consideration move the 
latter down until it comes into contact with either the 
substrate or another particle which lies immediately be- 
low it. At this point motion ceases. In other words we are 
looking at a model of BD plus surface relaxation where 
the latter is defined to be vertical movement downwards. 
At this stage the extent of the movement downwards is 
ill defined. It is obvious that such a model of BD plus 
vertical movement downward gives us exactly the same 
end result as (pure) RD. Hence the equation of motion 
for BD plus vertical movement downward (or pure RD) 
is Eq (|^). However, the interface width for such a situ- 
ation increases without limit and hence, after an infinite 



time, the width is infinite. Correspondingly the verti- 
cal movements downward must be infinite. What this 
shows us is that when we define the rules of movement 
for BD plus vertical relaxation, we have no idea of the 
extent of the vertical movement. However, the model, 
when evolved self-consistently according to these same 
rules, shows us that the motion in the vertical direction 
is unlimited. Furthermore, although vertical movement 
is occurring, the equation of motion, i.e. Eq (||), does not 
contain any differential operators, the structure itself is 
free of vacancies and the surface is not self-affine. 



1. Surface Relaxation Involving Vertical Motion Upwards 
Only 

Let us consider next a variant on this model but this 
time involving upward motion. In essence, the previous 
model started from (pure) BD allowed vertical motion 
downwards and ended up with (pure) RD. We will now 
consider the converse, i.e. we will start with (pure) RD 
allow vertical motion upwards and end up with (pure) 
BD. Thus the rules of motion are the following. Choose 
a random number and add one particle to the correspond- 
ing lattice point. If the added particle has no neighbour- 
ing points occupied by particles to a higher level, then 
leave it where it is. On the other hand, if the added par- 
ticle finds higher columns next to it, move it up the side of 
the highest such column until it is level with the top most 
point. Again at this stage we have no idea of the extent 
of the upward motion that will ultimately be involved in 
such a model. It is well known of course that (pure) BD, 
which is in essence the model we are describing, leads 
to a surface width which ultimately saturates. In turn 
this means that the upward motion is in fact finite and 
determined by the surface width. Furthermore the asso- 
ciated equation of motion for the surface is well known. 
However we deliberately refrain from writing it down at 
this stage for reasons that will be apparent later. Suffice 
it to say that our original starting point was (pure) RD 
for which the equation of motion is Eq (^). We then in- 
cluded what turns out to be finite upward motion. The 
latter can, apparently, be described by adding two dif- 
ferential operators to Eq (||) . One of these turns out to 
have the same form as that given in Eq (^. However the 
other cannot be written as a divergence term i.e. it is 
a non-conservative term and is in fact represented by a 
non-linear operator. 

Another interesting aspect of these two cases is the fol- 
lowing. At the outset it is not apparent in either model 
just what the extent of the vertical motion is. How- 
ever the models themselves evolve in such a way that the 
downward vertical motion is unlimited in extent whereas 
the upward vertical motion is limited. This shows that 
the rules of movement, which appear to be the 'mirror 
image' of one another, lead to models which evolve in 
such a way that up-down symmetry in the growth di- 
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rection is broken, i.e. motion downwards can occur to a 
significantly greater extent than motion upwards. 

2. Surface Relaxation Involving Limited Vertical Motion 
Downwards Only 

We now wish to consider a third model which, as we 
will show, turns out to be a hybrid of the previous two. 
When considering horizontal movement, we argued that 
the existence of universality classes, associated with sur- 
face relaxation processes, was in fact a consequence of the 
limited amount of horizontal movement that was permit- 
ted to occur. We will now examine the corresponding 
question of imposing a limited amount of vertical move- 
ment downwards. Thus we begin with a model which 
initially starts from (pure) BD and then allows a limited 
number of vertical moves downward (i.e. a maximum of 
100 lattice spacings say). Initially, when "growing" such 
a model on the computer, we will be dealing with essen- 
tially pure RD. This is because the interface width in the 
early stages of 'growth' will be much less than 100 lattice 
spacings. Hence initially we will have a SOS pure RD 
growth mode in which there are no vacancies present in 
the structure. For this stage of the growth the govern- 
ing equation will be Eq (^. However, as we know, the 
interface width grows without limit for (pure) RD and 
a stage will eventually be reached at which the interface 
width approaches and then exceeds 100 lattice spacings. 
Once this occurs a situation will develop where some of 
the downward vertical motion will leave vacancies in the 
structure, i.e. some of the downward motion will not 
be SOS. For such a situation we would expect from our 
earlier considerations that the form of the equation of 
motion would change and, by analogy with the previ- 
ous case, we might anticipate non-linear terms to appear 
in it. The hybrid nature of the present case from the 
earlier viewpoint of describing (pure) BD starting from 
(pure) RD is that only some of the particles having higher 
columns next to them are eventually moved up, and even 
then they are only moved part of the way up the column. 

The essential point here is the following. We started 
from a model with well defined rules of movement. In 
the early stages of growth, the resultant structure was 
defect free (i.e. no vacancies) and was described by a 
particular equation. However, as a self - consistent result 
of the rules of movement, the surface eventually evolves 
according to a different differential equation from the one 
describing the earlier stages of growth. During the sec- 
ond stage of growth defects (i.e. vacancies) do occur in 
the structure and the resulting equation would be antic- 
ipated to contain non-linear terms. It is important to 
stress that it is not the case that the equation of mo- 
tion is the same throughout all the growth, but with 
a cross-over from one regime of dominance to another 
regime. Rather different forms of equation are required 
to describe the evolution of the different stages of growth. 
What is more, if we increased the maximum number of 



allowed downwards vertical moves from 100 to say 1000, 
the time during which the first growth stage was gov- 
erned by pure RD would increase correspondingly. This 
shows clearly that the "cross-over" time between the two 
regimes is governed by the number of vertical moves that 
are permitted to occur. Given that (pure) RD is not 
a self-afhne surface, we can see that in general whether 
such a surface (i.e. self-affine) is ever seen to evolve in 
a finite time computer model can be crucially dependent 
on the nature and extent of the permitted movements. 

C. Heuristic rules 

As a prelude to describing the proposed rules we con- 
sider a simplified model which permits some insight into 
this question. Consider first pure random deposition in 
one dimension and the question of the change in the inter- 
face width following the deposition one monolayer. The 
expression for the final height (h) f reads 

{h)f = j[J2ih, + l + Vr)] = {hh + l (8) 

i 

where {h)i is the initial average height and L the 'length' 
of the substrate. 
Similarly 

('''h- = i[E('^' + 2/^^ + 1) + 2'7.(/i. + 1) + V^] 

i 

i.e. 

{lr?)f ^ {h\ + 2{hU + 1 + {r^f) (9) 

It follows from equations ^ and |^ that the final inteface 
width Wf is related to the initial inteface width Wi by the 
relation 

w}^{h')f-{h)}^wf + {rjf) (10) 

If we now imagine a relaxation process occuring in- 
volving the movement of particles downwards in which a 
fraction a of the columns lose a particle and a fraction a 
of them gain one, a similar analysis leads readily to the 
result 

wj = wf + (r/f ) +2a + 2a[(h,)g - (/i,)^] (11) 

where {hi) g{{hi) l) is the aveage height of the columns 
which gain (lose) a particle. It follows from Eq. [l^ that 
saturation will occur if the following relation is obeyed. 

{h^)L~{h.)g^^{vf) + l (12) 

Since the noise term (i.e. (r/f)) is fixed it is clear that it 
will be difficult, in general, to satisfy this equation on a 
layer by layer growth basis and that in general large fluc- 
tuations (or oscillations) about the 'equilibrum' width 
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will occur in computer simulations on finite size sub- 
strates. Another important feature that is apparent from 
examination of Eq. 1^ is the following. If, within a class 
of allowed moves, we permit only a fraction (3 of any given 
move to occur, we would replace a by the fraction a[3. 
This means that the entity on the left hand side would 
increase progressively as we decreased the value of (3. In 
turn this implies that the corresponding interface width 
would increase. Since allowed moves of a given class are 
presumed to be associated with a universality class of a 
particular type, this implies that the interface width of 
the latter can be increased continously simply by reduc- 
ing the probability of movement uniformly for all allowed 
moves in the class (clearly the range of possible values of 
/3, although undefined at this stage), can not be over the 
entire region < /? < 1 since e.g. = gives the model 
involving no relaxation i.e. exactly how small /3 can be is 
not known at this stage. Put alternatively, this suggests 
that interface width alone is a poor indicator of the uni- 
versality class since the prefactor of the term describing 
the time evolution of this width for such a class can be 
increased by simply varying the probability of occurence 
of all allowed moves in the class. It is revealing to pur- 
sue this argument further and to distinguish between the 
various types of move contained within a given class of 
moves. This can be done on the basis of the co-ordination 
number Vi of the particle prior to movement. For exam- 
ple in case of = 1 we could distinguish between the 
fraction of particles ai having a Ui = 1 and sat on top 
of a column, from the fraction having a Vi = 1 and 
sat next to a step edge. Yet again these are to be distin- 
guished from the fraction as of particles having a = 2. 
For generality we will assume that the fraction of parti- 
cles are ai up to a„. (In terms of the previous notation 
(a = ai -|- q;2 -t- . . . + The generalisation of Eq. |ll| 
can be readily shown to be 



w. 



+ i-nf) +2a + 2{ai[(/i,)g - 



- a2[{hi)2g - {hi) 



2L\ 



i/ng 



{h.U]} (13) 



where, in an obvious notation {hi) ng{{hi) ul) is the av- 
erage height of type n-columns that gain (lose) a parti- 
cle. Once again the effect of reducing all the ai by the 
same factor /3 will give the immediate generalisation of 
the result described earlier. However it is clear from the 
present equation thath the functional form of the equa- 
tion remains unchanged only when all the ai are reduced 
in this way. As opposed to this, if the ai are all reduced 
by different factors j3i, then the functional form of the 
equation is changed. For such a situation it is far from 
obvious whether the resultant class of moves belongs to 
the same universality class as for [3 — 1 (or for a uniform 
value of /3 ). In the extreme case that /3i = for some 
of the moves, but exists for others, we would clearly an- 
ticipate a different universality class. However, even for 
the latter situation, it is not clear whether increasing the 



j3i values which were previously identically zero to ex- 
tremely small values would change the universality class 
type, i.e. the stability of the latter to small changes of 
movement type is not presently known. These questions 
will be investigated in subsequent publications. 



D. Heuristic arguments for the differential operators 
appearing in the stochastic differential equation 

To date, by considering simple "growth" mechanisms, 
we have argued that the question of the nature of the 
evolving surface, and the precise form of the stochastic 
equation describing its evolution, is determined essen- 
tially by the nature and extent of the 'horizontal' and 
'vertical' movements that are permitted to occur. In the 
present section, based on our considerations in the pre- 
vious one, we will present simple heuristic arguments for 
the form of the differential operators to be associated 
with such movements. 

Before proceeding to detail it is informative to consider 
the range of possibilities that could occur for an evolv- 
ing surface. Clearly the associated interface width could 
either saturate or not. For the latter, there are then 
two possibilities - namely that the associated stochastic 
growth equation contains differential operators or it does 
not. If on the other hand the interface width does satu- 
rate it is almost certain that the associated growth equa- 
tion does contain differential operators. The question is 
which ones? 

We form next a set of rules for which differential oper- 
ators are to be associated with which situation (see also 
Figure]^ ). In the latter, following conventional wisdom, 
we have considered four types of differential operators 



only namely (f )2 ^(£)2 and |^ (As should 

be clear from Fig ^ the effect of such operators on a 'sym- 
metric hill' is symmetric - hence as long as the prescribed 
rules of deposition and surface relaxation possess this left 
- right symmetry one would expect such even order op- 
erators). 

Our suggested rules (based simply on heuristic argu- 
ments and arising from a consideration of the different 
types of interfaces that evolve from simple growth mod- 
els) are the following. 

For surfaces whose interface width saturates; 

• Movement vertically downwards (further refine- 
ment of these rules are presented later) 



1 



employ the operator if the rules of 

movement include all possible downward mo- 
tion 

& 



employ the operator — if the rules 
movement exclude certain Vi — 2 particles. 



of 



Movement vertically upwards (further refinement 
of these rules are presented later) 
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1. employ the operator +7(f|)^ if the rules of 
movement involve the creation of vacancies 
in the structure and correspond to maximal 
movement (with 100 % probability) of parti- 
cles up the side of columns but do not involve 
movement to the top of the adjacent column 
(i.e. non SOS models). 

2. employ the operator ^-§^{^Y if f^'^ rules of 
movement do not involve the creation of va- 
cancies in the structure and correspond to the 
movement of particles (with 100 % probabil- 
ity) upward from the top of one column to the 
top of an adjacent column (i.e. SOS models). 
Note that upwards motion can, on its own, 
destabilize a surface, and for SOS relaxation 
needs to be conterbalanced by an equivalent 
downwards motion 'current' 

For surfaces whose interface width does not saturate 
there are at least two distinct possibilities. Either the 
growth equation does not contain differential operators 
(e.g. pure RD or, as we will argue in a later paper, pure 
shading) or it contains an instability. To see how the lat- 
ter could arise consider the case of pure RD plus (limited) 
surface relajxation (i.e. the deposition process modelled 
by Family Bl) As described earlier, the evolution of this 
surface in d' = 1 is described by Eq (^ and snapshots 
of this surface at successive stages of growth is shown in 
figure 1^, i.e we have a smooth surface, in this instance, 
in which the rules of movement lead, self-consistently, to 
a succession of small horizontal and small vertical down- 
ward motions. 

Suppose we have a video film of the evolution of this 
surface as it runs forward in time and (unknown to the 
viewer) we instantly stop the forward running of this film 
and instead run it backwards. Mathematically the re- 
sulting unwinding of the surface relaxation process cor- 
responds to time reversal and is equivalent to replacing 
t by —t in Eq (^. One consequence of the latter is that 
we would replace u hy —u va Eq (^. Clearly this would 
now correspond to small vertical upward motion (see fig- 
ure H for a graphic illustration of this) Eventually such 
upward motion would 'undo' the effects of the previous 
surface relaxation processes and we would in effect move 
back towards pure RD - i.e. the surface would become 
unstable as far as saturating the interface width was con- 
cerned. In other words, the existence of a term — A^j^ 
could destabilize the surface growth. (This result can 
also be readily established via the formal machinery of 
the group renormalisation approach.) It should be clear 
from these considerations that reversing the sign of the 
four differential operators we discussed earlier could also 
result in instabilities in the associated differential equa- 
tion. Hence one has to examine carefully the permitted 
movements in order to ensure that such movements do 
not ultimately lead, via self-consistency, to terms in the 
associated differential equation which make its solution 



unstable. An example of this is given later. 

E. Comparison Of The Heuristic Rules With 
Computer Models Reported In The Literature 

Most computer models reported in the literature per- 
mit the deposited particle one move in the horizontal 
direction. Which direction the particle moves vertically 
is also determined by the rules of movement. In certain 
cases both upward and downward motion is permitted, 
whereas in other cases upward (or downward) motion 
only is permitted. For many cases reported in the liter- 
ature the rules governing surface relaxation (which are 
almost invariably of the SOS type) are not even stated 
in this simple form, rather, in an effort to mimick the 
physics, the rules are defined in terms of the coordina- 
tion number Vi of the particle to be moved. In particlular, 
when Vi > 2 movement is not usually permitted unless 
the coordination number increases - i.e. in those situa- 
tions where the coordination would be the same at the 
end of the movement as it was prior to movement (the 
so-called 'tie' situation) no movement is permitted. On 
the other hand where the coordination number increases 
via the movement, the rule could either be 

(a) move to the nearest neighbour site that increases 
the value of Vi 

or (b) move to the nearest neighbour site which gives the 
maximum value of Vi. 

The subtlety of these rules is that, because of the topol- 
ogy of the surface, (a) and (b) define different types of al- 
lowedimoves in different dimensions. For example Kotrla 
et al.Q have shown that both (a) and (b) lead to vertical 
motion downwards in c?' = 1. However in d' = 2 rule (a) 
leads to predominantly downward motion whereas rule 
(b) leads to both downward and appreciable upward mo- 
tion. (Further clarification of this point in quantitative 
terms will be given in a subsequent paper.) 

We will now consider various models from this viewg 
point. For the well known d! = 1 case studied by Familytl 
the 'tie' situation occurs frequently for = 2 and the 
rules of movement permit the particle to move down only 
under these circumstances. Correspondingly the result- 
ing interface is smooth - hence we are led to the conclu- 
sion that if all moves downward are permitted (including 
the tie situation) the resulting interface is smooth and 
the associated relaxation downward is described by the 
operator v^^- As opposed to this it is found that if, 
in d! = 1, moves downwards for the tie situation when 
Vi — 2 are forbidden (Figures I and I) and downward mo- 
tion only allowed if the coordination number is increased, 
we obtain a rough interface (Fig M whose properties are 
described by the operator — k^^^B The generalisation 
of these same rules to c?' = 2 presents an interesting sit- 
uation in that for case (a) discussed above, this leads to 
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predominantly downward motion only and is hence de- 
scribed by the operator — K^^-Eli However for case (b), 
significant upward motion is also allowed which is essen- 
tially SOS in nature. The coresponding interface (which 
is rough) will be described by an equation in which the 

downward motion is represented by the operator —n 



whereas the upward motion is decribed by A^3-(^)^. 
This is in-agreement with what is found in the computer 
modellingp'LJ. Hence we here have a situation where the 
equation of motion, because of the permitted rules of 
movement, is governed by the following equation 



dt 



in d' = 1, but is described by the equation 



dh 
dt 



^ 9^ ^dh^2 ^ ^ 
dx^ dx 



(14) 



(15) 



in d! = 2. We would argue, in contrast with statements 
in the literatur^ that this is not an unusual situation in 
physical terms but is merely a consequence of the fact 
that the permitted movements in dl = 1 and d! — 2 are 
different. Equally we would argueQ that it is not true 
that the equation of motion is the same in d' = 1 and 
d! = 2 but that the crossover has not been seen iri-|d' = 1 . 

Consider next the model of Lai and Das SarmaS These 
authors considered |a similar model in d' = 1 to the onje 
of Wolf and VillaintI and Das Sarma and Tamborenea.tl 
However in the case of tie with Vi = 2 the particle was 
moved to the nearest neighbour site with the smaller 
height difference (i.e. sometimes upwards and sometimes 
downward motion occurred). For the case of a tie with 
Vi = 1 the particle was allowed to diffuse to higher bond- 
ing within a distance I - i.e. for d' = 1 downward mo- 
tion could also occur for this situation. Given that these 
rules are all SOS type and that significant ammounts of 
upward and downward motion is occurring, we would ex- 
pect that the resulting equation is the fourth order non- 
linear equation appeariiLg in Eq ( p^ . This is precisely 
what Lai and das SarmaQ found. 

In view of our earlier comments we would envisage an 
interesting situation developing if, in the case of tie with 
Vi = 2, the downward motion was forbidden and only 
the upward motion was allowed. For such a situation 
we could anticipate an instability could develop in the 
system since such upwards motion is not counterbalanced 
by an equivalent downwards motion 'current'. [This is 
precisely what happens in the model of Park et al.EI where 
no saturation of the interface width was found for systems 
larger than a critical size. 

As a final example we consider model 1 in d' = 2 of 
Kotrla et. alD. This model is a rStraightforward gener- 
alisation of the Wolf and VillainO and Das Sarma and 
Tamboreneafl model in that if there is a neighbouring 
site with a higher value of Vi the particle relaxes to it 
in either the upward or the downward direction In the 



case of a tie the particle remains where it is. It should 
be obvious from our previous discussion that once again 
we have SOS relaxation in the upward and the downward 
direction and that the fourth order non linear equation 
i.e. Eq ( p^ ) will describe the evjolution of the surface. 
This is exactly what Kotrla et al.cl found. 

Given the success of our heuristic rules we could then 
invert the problem of predicting the permitted move- 
ments to be associated with a given stochastic differential 
equation. Consider for example the equation 



dh 
'dt 



F 



d^h 
dx^ 



i9h,2 



(16) 



The terms F + imply that the initial deposition is 
SOS and is pure RD. The second term represents move- 
ments horizontally and vertically downwards. Similarly 
the term (fj)^ implies relaxation vertically upwards with 
the particle being moved up the side of a column until it 
was (at most) placed level with the top of it (i.e. non SOS 
vertical motion which creates vacancies in the system). 
In summary, we would expect that Eq (^6|) represented 
RD following by surface relaxation in which the process 
involves both downward and upward motion. It comes 
as somewhat of a surprise therefore to find that Eq ( p^ ) 
is apparently applicable in d' = 1, to pure BD. If any- 
thing we might have expected that the latter corresponds 
to (pure) RD followed by limited vertical movement up- 
wards along the side of higher (neighbouring) columns - 
i.e. one might have expected (pure) BD to have obeyed 
the equation 

dh ,dh.o 

= ^ + 

A possible resolution of this paradox would be (as in- 
dicated by our previous reasoning) that (pure) RD fol- 
lowed by surface relaxation in both the downward and 
upward direction is mathematically indistiguishable from 
pure BD. This problem will be addressed in future publi- 
cations. For the present we note that, in graphical terms, 
what we are suggesting is that the situation depicted in 
Fig (a) (appropriate to 'pure' BD) is ultimately physi- 
cally indistinguishable from that depicted in Fig ^ (b) i.e. 
RD plus surface relaxation in the downward and upward 
direction. Alternatively we could argue that we need 
both a downward and an upward 'current' of particles to 
obtain a stable interface. 



III. CONCLUSION 

In the present paper we have presented arguments to 
the effect that observation of universality classes in com- 
puter simulations of surface growth is a consequence of 
the limited number of horizontal moves that a given par- 
ticle is permitted to make. (Conversely, if the particle 
was permitted to make any number of moves up to a max- 
imum value determined to be equal to the substrate size. 
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then interface roughness would never develop in the sense 
that growth would be of the layer by layer type). Given 
the limited number of horizontal moves, the universal- 
ity class is then determined by the nature and extent of 
the vertical motion. This ultimately is controlled by the 
rules associated with the coordination number Vi taken 
in conjunction with the dimensionality d' of the surface. 
In particular rules may be defined which in effect allow 
for example 

(a) downward motion only in both d' = 1 and d' = 2 

or (b) downward motion in d' = 1 but both downward 

and upward motion in d' = 2 

If downward motion only is allowed and this is of the 
SOS type (i.e. a particle has another particle beneath it 
at the beginning and end of the motion) then the degree 
of the differential operator i.e. or is determined 
by the rules of motion associated with the tie situation. If 
motion downwards is allowed in the tie situation, then the 
second order differential operator is appropriate since this 
corresponds to all moves downwards being permitted. On 
the other hand, if motion downwards is forbidden in the 
tie situation, then the fourth order differential operator 
is applicable since such an operator is appropriate for 
the situation where only a fraction of the moves that can 
occur are permitted to occur. Similarly if both upward 
and downward motion is allowed in the tie situation the 
fourth order equation is again applicable. 

In the case of upward motion the form of the differ- 
ential operator will again be determined by the rules of 
movement. If the upward motion is of the SOS type then 
it is decribed by the operator ^(f|)^- However if it is 
non SOS type (i.e. it leads to the creation of vacancies) 
it will be described by the operator (f^)^. 

We have shown, in the body of the text, that these 
heuristic rules are consistent with many standard results 
reported in the literature. Furthermore we have pre- 
sented arguments to show that the nature and extent 
of the permitted vertical motion can lead to situations 
where, for example, initially the growth contains no va- 
cancies and is described by a particular differential equa- 
tion. However, eventually vacancies will enter the struc- 
ture (as a result of the rules of movement) and subse- 
quently for this situation, the growth is described by a 
different differential equation. 

Further substantiation for the validity of these argu- 
ments will be presented in a series of subsequent papers 
where we will also address the question of the apparent 
anomaly of the equation of motion for pure BD, and the 
question of the universality class when the /3j is chosen 
to have different values for the different move types i. 
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FIG. 1. DifFcrcntial Operators acting on the function h{x) 
increases the particles at the side of the 'hill' but does 
not 'transport' them to the top of it. 

removes particles form the top of the hill and redis- ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^H 

tributes them to the sides and bottom of it. 

- 1^ is similar to the previous case. FIG. 3. Surface of the Wolf and Villain Modell (taken from 

+-§^{^f moves particles from the base of the hill to the Ref. 1) 
top of it. 
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FIG. 4. Growth rules of Family 




FIG. 5. Growth rules of Das Sarma and Tamborenea 




FIG. 6. Growth rules of Wolf 




FIG. 7. Different Rules With The Same Result 
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